ON LIOUVILLE SYSTEMS AT CRITICAL PARAMETERS, PART 1: 

ONE BUBBLE 



CHANG-SHOU LIN AND LEI ZHANG 

Abstract. In this paper we consider bubbling solutions to the general Liou- 
ville system: 

, A J hie"' \ 
(0.1) AgUi + Y,aijPj \— — F-1 =0 in M, L...,n(n>2) 
j=l \Ihje"j J 

where {M,g) is a Riemann surface, and A = {aij)nxn is a constant non-negative 
matrix and pj — > pj as A; — 5- <». Among other things we prove the following sharp 
estimates. 

(1) The location of the blowup point. 

(2) The convergence rate of Pj — pj, j = l,..,n. 

These results are of fundamental importance for constructing bubbling solutions. 
It is interesting to compare the difference between the general Liouville system 
and the SU (3) Toda system on estimates (1) and (2). 



1. Introduction 

Let {M,g) be a compact Riemann surface whose volume is normalized to be 1, 
hi,...,h„ be positive functions on M, pi, ..,p„ be nonnegative constants. In this 
article we continue our study of the following Liouville system defined on (M, g) : 

(1.1) ^^Ui+f^pJaiJ{-^^^^—--\)=0, /G/:={1,..,4 

where dVg is the volume form, A = is a non-negative constant matrix, is 
the Laplace-Beltrami operator (— > 0). When n = I and an = 1, equation (11.11) 
is the mean field equation of the Liouville type: 

Therefore, the Liouville system (11. Il l is a natural extension of the classical Liou- 
ville equation, which has been extensively studied for the past three decades. Both 
the Liouville equation and the Liouville system are related to various fields of ge- 
ometry. Physics, Chemistry and Ecology. For example in conformal geometry. 
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when p = Stt and M is the sphere S^, equation ( I1.2I) is equivalent to the famous 
Nirenberg problem. For a bounded domain in M? and « = 1, a variant of (11.21 ) can 
be derived from the mean field limit of Euler flows or spherical Onsager vortex the- 
ory, as studied by Caglioti, Lions, Marchioro and Pulvirenti f6','7'|, Kiessling |28], 
Chanillo and Kiessling [9J and Lin |34]. In classical gauge field theory, equation 
(II. lb is closely related to the Chem-Simons-Higgs equation for the abelian case, 
see [51 |24l |25l im. Various Liouville systems are also used to describe models 
in the theory of self-gravitating systems [Yj, Chemotaxis lITTl l27l . in the physics 
of charged particle beams ||4]|20l|29]|30l, in the non-abelian Chem-Simons-Higgs 
theory I2T1 l26l l45l and other gauge field models Il22]|23ll3n . For recent develop- 
ments of these subjects or related Liouville systems in more general settings, we 
refer the readers to ||2] S [E] |T2 112 [B] III] [11] |32] |32 134] |^ 
|43]|Slll6l|43 and the references therein. 

For any solution u of (11.21) . clearly adding any constant to u gives another solu- 
tion. So it is nature to assume u ^H^{M), where 

H\M) = |m G L^{M) |VgM[ G L^{M) and J udVg = o| . 

Corresponding to (11.11 ) we set 

H^" =h\m) x---xH\M) 

to be the space for solutions. For any p = (pi, • • • ,p„), p; > 0(/ G / = {1, ...,«}), 
let <Pp be a nonlinear functional defined in H^" by 

*p(«) = \L '''^ I VgMrVgM,t/Vg-£pylog / hje'^idYg 
^ijei ■J'^ jei ■J'^ 

where (a'^ )„x« is the inverse of A = {aij)nxn- It is easy to see that equation (11.11 ) is 
the Euler-Lagrangian equation of <l>p . 

If the matrix A satisfies the following two assumptions: 

{HI) : A is symmetric, nonnegative, irreducible and invertible. 
{H2) : a" < 0, V/ G /, a'^ > 0, V/ / j, £ a'j > 0, V/ G /, 

the authors prove in |I3F| that for p satisfying 

(1.3) SttA^^P; < £ aijPiPj < Sn{N+l)Y,Pi, 

iel ijel iel 

there is a priori estimate for all solutions u to (11.11 ). and the Leray-Schauder degree 
dp for equation (11.11 ) is 

dp = ^ [i-XM + 1)...{-Xm+N)^ if (O holds 

where Xm is the Euler characteristic of M. Moreover, if p*^ tends to the hyper- 
surface {p; ^TlNY^ieiPi = Hi.jei'^ijPiPj}^ there exist exactly disjoint blowup 
points (see l,38j ). 
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The proof of the a priori bound in 11381 reUes on the sharp estimate for a sequence 
of bubbling solutions to (ll.ll ). Let be the blowup solutions corresponding to p*^ 
and B{pt,8{)) (t = l,..,N) be disjoint balls around distinct blowup points in M. 
Then under assumptions (HI) and {H2), the behavior of u'^ around any p, is fully 
bubbling, that is, the maximum values of any components of u'^ in any of the balls 
are of the same magnitude: 

max = max + 0(1), V/, /' G /. 

B{p„So) B{p„So) ' 

Moreover, after a suitable scaling around each blowup point pt, converges to an 
entire solution f7 = (L^i, ..,^7„) of the following Liouville system: 



(1.4) 



^Ui + l:]^^aije^J =0, in 



e^' < oo Ui is radial, V/ G /. 



One may expect the limiting entire solution to be different around each blowup 
point, however the authors proved that U is independent of blowup points, and only 
depends on the ratio of — pi : p| — p2 : ..p^ — p„ (see |37 1). Naturally it leads to 
the question: how to construct bubbling solutions with the help of this information? 

In this paper and subsequent ones, we are devoted to study the bubbling phe- 
nomenon of Liouville systems: how to accurately estimate the bubbling solutions 
of dl.ll ) and how to construct them. These are quite challenging analytic problems. 
In general, blowup analysis for a system of equations is much harder than that for 
the single equation. One reason is that the Pohozaev identity, a balancing condition, 
is no longer so powerful as in the scalar case. Another reason is that there are too 

many entire solutions: the parameter a = (ai , a„) (a,- = — / e"'), which rep- 

resents the energy of the entire solution, forms a submanifold of n — \ dimension. 
However, for the Liouville equation, the energy is just one number: J^2 e" = Sn. 

In this article we consider the case of one blowup point, and always assume 
(HI) only. Let p = (pi, ...,p„) and 

(1.5) A7(p) = Stt^p,- - £ aijPiPj 

ieJ i.jeJ 

for any J Cl := {l,...,n}. Define 

n = {p; A/(p) = and Ay(p) > for all C / C / }. 

Note that if {H2) also holds also, then A/(p) = implies A7(p) > for all 7 C / ( 
see 1.38,1 ). For any p we define (nii, ...,m„) by 

(1.6) = ^Y^c^ijPj- 

The quantity m, can be interpreted by the entire solution U of (11.41 ). In fact 

Pi = / e^'dx, i G /, 
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and 

(1.7) Ui{x) = —m, \og \x\ + 0(1), for \x\ near infinity. 

The integrability of e^' implies m, > 2 for all /. On the other hand A/(p) =0 can 
be written as Y^iei ("^! — 4)p; = 0. Thus either min{mi , . . . , m„ } < 4 or m,- = 4 for all 
/ G /. We also note that (11.71 ) imphes that Fi is a smooth submanifold because the 
normal vector at p G Fi : 

(L "UPj - 471, £ a„jpj - 471), 
jei jei 

has all its components positive. 

The asymptotic behavior of Ui{x) shows that the decay rate of e^'^^^ is 0{\x\^'"), 
where 

m = min{mi, ...,m„}. 
In this article we define 2 G Fi with m = 4, i.e. m, = 4 for all /. Thus the decay 
rate of e^' for p = 2 is 0(|xj^^). The difference on the decay rate for Q and p ^ Q 
will have great effects on bubbling analysis later. 

Let = (mj,...,m^) be a sequence of blow up solutions to (11.11 ) with p = p*^ 
such that p*^ — > p G Fi. The point Q defined above is of particular importance, 
the readers will see that in our main theorems below, the asymptotic behavior of 
blowup solutions, the nature of A/(p*^) and the location of blowup point are all 
significantly different depending on p = 2 or not. 

Let p be the blowup point of and 

(1.8) Mk = max f ^^(x) -log / he^'dvA, 

B{p,S) V JM V 

(1.9) £^ = f>-I«*. 

Since p*^ — )• Fi, there is only one blowup point p. It is easy to see that fully 
blows up at p (see Lemma [6?!] ): 

(1.10) max (4W-log / V'JK) =Mi + 0(l), / G /. 

B{p,So) JM 

Our first result is on the location of the blowup point p. Let pt ^ phe where 
the maximum of {u\ , . . , } is attained, then we have 

Theorem 1.1. Let p*^ ^ p G Fi and all pf — Pi have the same sign. 

(1) ifp / a then 

(1.11) £ (v{\oghi){pt)+2KnN,Y{pk,pk)^Pi = 0(£f-'). 

(2) Ifp = Q, then 

(1.12) (v(loghi){pk) + SKViY{pt,Pk)^Pi = 0{ei\oge,'). 

where Vi denotes the derivative with respect to the first variable, and 
7(x,3') stands for the regular part of the Green's function. 
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Our second result is about the decay rate of A/(p*^). To state the result, we need 
to define the following quantity D, (/ G / = {1, ...,«}): 

(1.13) Di=\im(5t'"-'^ I M4g2™(G{^-,p)-r(p,p))^v) 

So^O^ " 2k Jm\B{p.So) hi{p) 

provided that m < 4. The limit is well defined if m < 4, see section 6. 

Theorem 1.2. Suppose p*^ ^ p € Fi and p ^ Q, if all — Pi have the same sign, 
then 

Az(p*^) = 8;^2£(e^'D,■ + o(l))£^^ 

ieh 

where h is a subset of I where m, = mfor all i h, c,- is a constant determined in 
KM . o{l) ^ as k ^ oo. 

If M is a flat torus with fundamental cell domain H C M^, then D, can be written 

as 

Di = ( / —dx — lim / ' ^ ''^} dx 

2k \Jm.i\Q.\x- p\"' Sfi^0jQ.\Bii^^\x- pf 

where 

HAx.p) = ^2™(r(.,p)-r(P,P)) _ / g /. 
hi{p) 

See O and lITSl for related discussions. 

Remark 1.1. The assumption that all pf — p,- have the same sign seems due to 
some technical difficulties. When n = 2 and p ^ Q, this assumption is not needed 
for both Theorem \l.l\ and Theorem \1.2\ 

Theorem 1.3. Suppose p* — > p G Fi and p = Q. If all pf — Qi have the same sign, 
then 

(1.14) Ai{p^) = -l6Ti'{Y.^ie'< + o{l))el\oge^\ 

iei 

where 

1 



^ ^ A \oghi{p) - 2K{p) + 871 + |V \oghi{p) + %nViY{p,p)\ 

Ci is determined in H6. 14i . 

Important information on bubbling solutions can be observed on the two cases: 
p i^Q and p = Q. Theorem 11.21 which is on p 7^ Q, has its leading term in A/(p^) 
involved with global information of the manifold, while the leading term in Theo- 
rem [T31 which corresponds to p = 2, only depends on the geometric information 
at the blowup point. Moreover, the error terms in Theorem 11.21 and Theorem 11.31 
respectively also indicate the different asymptotic behaviors of blowup solutions 
near the singularity. All these differences in the two cases will lead to separate 
strategies in the construction of bubbling solutions in forthcoming works. 

Since Liouville systems and Toda systems share a lot of common features, it 
is informative to compare our main theorems with the ones for the SU (3) Toda 
system. First, the location of the blowup point in Theorem II. H is a critical point of 
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a combination of log/j,, P; and 7 (Vi 7 vanishes if the Riemann surface has constant 
curvature). However for the SU (3) Toda system, the blowup point p is a. critical 
point of both log/zi and log/z2, i.e. p satisfies (see |[36l ) 

Vh,{p) = Vh2{p)=0. 

Second, for the SU (3) Toda system, the convergence rate of — p, is estimated to 
be 

where hi is the term in (I1.141 i. Nevertheless our result in (I1.14t is again a com- 
bination of the biS. The comparison of the results reflects some major difference 
between the Toda system and our Liouville system: 

(1) The dimension of kernel space of the linearized operator at an entire solu- 
tion is 8 for SU (3) Toda system, and is 3 for our Liouville system. 

(2) Thesetr:= {(pi,...,p„); pi = J^2e"', (mi , is an entire solution } 
is only a point for SU (3) Toda system, while for the Liouville system it is 
a (« — 1 ) dimensional manifold. 

As far as the blowup analysis is concerned, our Liouville system has disadvan- 
tages in both respects, as the kernel space is too small and F is too large. For a 
sequence of bubbling solutions, it is extremely difficult to pin-point suitable ap- 
proximating solutions from F, because at the beginning, the local energy of bub- 
bling solutions could be estimated in some rough way. This rough estimate of 
the local energy leads to a small perturbation of global bubbling solutions. This 
perturbation on global solutions, albeit small, has a non-negligible effect on the 
approximation of blowup solutions. This difficulty is particularly evident when 
we study bubbling solutions with multiple blowup points in [39|. Therefore our 
method to obtain those sharp results is different from ones in Chen-Lin [ 12] for 
the mean field equations and Lin- Wei-Zhao |[36l for the SU (3) Toda system (The 
methods in IIT2I and ||36]| are similar). 

The organization of the paper is as follows. In section two we first prove a 
uniqueness theorem for globally defined linearized Liouville systems. This re- 
sult plays a central role for the delicate blowup analysis in sections three to five. 
The main idea of the proof uses a monotonicity property of solutions and we in- 
troduce a way to use maximum principles suitable for Liouville systems. In the 
second part of section two, we study the asymptotic behavior of global solutions to 
the Liouville system on and obtain some Pohozaev identities. In section three 
and section four we obtain a sharp expansion result for blowup solutions around a 
blowup point. Then in section five, for local equations we use Pohozaev identity 
to determine the locations of blowup points. Then in section six we return to the 
equation on manifold and compute the leading term for p*^ ^ p in both situations 
ip ^ Q or p = Q) and complete the proofs of the main theorems. 

Acknowledgement Part of the paper was finished when the second author was 
visiting Taida Institute for Mathematical Sciences (TIMS) in March 201 1 and June 
2012. He is very grateful to TIMS for their warm hospitality. He also would 
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2. Preliminary results on the Liouville systems 

In this section we prove two theorems on the Liouville systems with the matrix 
A satisfying {HI). They are important for the blowup analysis and the computation 
of the leading terms of p — >• Fi in this paper and p — )■ F^r in the forthcoming work 
11391. 



2.1. A uniqueness theorem for the linearized system of n equations. In the first 
subsection we prove a uniqueness theorem for the linearized system of n equations. 

Theorem 2.1. Let A be a matrix that satisfies (HI), u = {ui,..,Un) be a radial 
solution of 

-Aui = 1 aije"' inR^, iel = {\,..,n} 
kie"' <oo. 

Suppose (j) = (0i,..,0„) satisfies 

n 

-A(l)i= Y,aije"J^j, iel. 

(1) 

(2.1) |0,-(x)| <C(l + |x|)^ xGM^, 

for some T S (0, 1) and 

0,.(O) =0, / G /. 
Then there exist ci,C2 € M such that 

0,-(x) = ciM-(r)— +C2M-(r)— in M?, i^l. 

(2) If\^i{x)\ <Cfor all X gM}, then there exist cq,ci,C2 £R such that 

0,.W=CoK(r)+2)+c,.:(r)^+C2«;(r)^, M^ iel. 

(3) If (pi{x) = 0{\x\'^) nearO and\(pi{x) \ <C{l + \x\)^^'^' for some Sq > 0, then 
^i = 0. 

Before the proof of Theorem 12. II we first establish 
Lemma 2.1. Let A and u be as in Theorem \2.1\ let = (Oi ,..,<!>„) solve 
^'/{r) + ^^'i{r) - ^<I>,(r) a,ye"^<I>y =0, < r < oo, 

|<J>/WI <Cr/(l+r)^ for some £0^(0, 1), V/ € /. 
Then there exists a constant C such that ^i{r) = Cu'i{r)for i € /. 



(2.2) 
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Proof of Lemma HiT} 

The proof is in two steps. First we show that under the assumption of 0, at 
infinity we have the following sharper decay estimate: 

(2.3) |<I);(r)| <Cr(l+r)-^ < r < oo, / g /. 

Indeed, let V/,(r) = <I>,(r)/r. By direct computation we see that ijf = , \j/„) 
satisfies 

(2.4) \i//ir) + ^ V>/(r) + £^,7^"^ Vi^; = 0, < r < 00. 

Clearly in order to show (I2.31 l we only need to show | V/, (r) | < Cr^^ for r > 1 under 
the assumption that < Cr^^' for / G / and r > 1. Let 

{t ) = \j/i{e') and m,- {t) = Ui{e*) + 2t, 

it is easy to see that = (1^1 (f), .., satisfies 

(2.5) i?>/'(0 + 2i?>/(0 + £a,7e'^Wv/y(0 = 0, — < ? < 00 

ye/ 

and our goal is to show 

(2.6) Ut) = 0{e-^') 
knowing v>,(f) = 0{e^^') for t >0,i £l. Set 

h = ^T,aij[ f-"', / = min{/i,. .,/„}. 

2n Jm.2 

By Lemma [2!2] below / > 2. Let 

ye/ 

Then ^ ^ 

V>',-(0=Co + Cif'-2' + - / hi{s)ds - -e-^' I e^'hi{s)ds. 

2 Jo 2 Jo 

Using the asymptotic rate of hi{t) at infinity we further have 

Yi{t) = {Co + - / hi{s)ds)+C,e-^' + 0{e^-^-'-'°^'). 
2 Jo 

Since Yi{t) tends to as f goes to infinity we know \j/i{t) = 0{e^^') if / + £0 > 4, in 
which case (12.61 ) is established. Otherwise we obtain Yiit) = 0{e^'^^^^^^'). In the 
latter case, we apply the same procedure to obtain a better decaying rate of ) at 
infinity. After finite steps, (I2.61 l is established. 

In the second step we complete the proof of the Lemma [TT] By way of con- 
tradiction we suppose there is a solution <I> = (<I>i, ..,<!>„) that satisfies (12.21 ) and 
<I> is not a multiple of / = {u\{r), ..,u'„{r)). Let i//^ = —u\{r)/r, then clearly both 
1//* = ()//{',.., 1//^) and Y = satisfy (12.41 ). We verify by direct computa- 

tion that 

f e'''\lf%o)o^ds=- f {e''')'o'^do = -e'''^''V + l f e^'iodoO. 
Jo Jo Jo 
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Note that the last inequahty is justified by Mj(r) < for r > and / G /. Also, since 

poo 

e"' < Cr^^^^ for some 5 > and r > 1, / e"'\lf^(a)o^do < oo. Based on the 

Jo 

computation above we set 

S= = + v/,)and ^ e"'i///(a)aVa > 0, Vr > 0, 

Let Y = {yi,.., Yn) and = {Yi, -jYn)- We first observe that V'f (0) > 0- Sup- 
pose Y ^ we can assume that Yi{0) / and |i/>,(0)| < i//f (0) for all / G /. 

From the definition of S we immediately see that € 5. Moreover, since | Yi{r) \ < 
Cr^ near infinity, we can choose |?| < 5 with 5 small so that all |f| < 5 belong to 
S. Another immediate observation is that S has a lower bound. Indeed, for T suf- 
ficiently negative, Yii'^) < 0, which is impossible to have /q e"' i/^f (a)a^<ia > 
for r small. 

Let r = inf5 and let f„, — )• r+. The sequence y/^"' obviously converges to a 
function Yh which is just Y? + T{Yf - ¥= iWi^ -lYn) satisfies the property 

(2.7) [' e"-^''^Yi{(^)(^^da>0, Vr > 0, i = \,..,n. 

Jo 

On the other hand, from the behavior of xj/ and at infinity (both are 0{r^^)) we 
immediately observe that 



/ e'"^''\i{r)r' dr <oo. 
Jo 



In regard to (I2.41 i we have 

(2.8) r^V^' W = - / £a,7e"' V^yC^^^^c^> < r < oo. 

JO j 

From ( 12.71 ) and (12.81 ) we see that Yi non-increasing. Since we have known that 
IV^K'')! < Cr^^ near infinity we have 

(2.9) '/^,('')>0. Vr>0, /G/. 

It is not possible to have all i//, (0) = because this implies Yi = 0, a contradic- 
tion to the assumption that Y is not a multiple of y^- Therefore without loss of 
generality we assume Yi (0) > 0- Then we further claim that Yi is strictly decreas- 
ing for all / G /. Indeed, let /i = {7 G /| fliy > 0}, for each j G I\ we use (12.91) and 
Y\ (0) > to obtain 

r^Y'iir) < - ( ^i/^'"' Yi<^^d<y < 0, < r < 00. 

JO 

Therefore for each j € h, Yj is strictly decreasing, which immediately implies that 
i//y(0) > 0. We can further define I2 = {i G /| atj > for some j G /i .}. Then the 
same argument shows that Yi is strictly decreasing for each / G h as well. Since 
the matrix A = {aij)nxn is irreducible, this process exhausts all / G /. 

( 12.81 ) yields Yii^) ^ —Cr^^ for r > 1 and / G /. Then by using hm,.^c» Yii^) = 
we further have 

(2.10) Yi{r)>Cr-^, r>l. 
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Then it is easy to see that for t = T — e with £ > small, we also have 

r e"'\l/l{a)a^da>0, forallr>0, 
Jo 

a contradiction to the definition of T. Lemma |2?T] is proved. □ 

Proof of Theorem l2.lt We first prove the third statement. The following function 
plays an important role: Let / = (/i , ..,/„) = {u[,..,u[^). Then 

(2. 1 1) - Afi = £ aije"^fj - ^fi, i=l,..,n. 

Let 0*^ = .., 0,^) be defined as 
1 f^^ 

^i{^) = — / ^i{r CO?, d,r&ind)co&kddd, iel. 
Ik Jo 

Then 0*^ satisfies 

(2.12) _A0f = £a,-/'^0}-^(/)f, iel, k = 2,.. 

j 

Clearly 0f (r) = o(r) near and ^^{r) = 0{r^^^) at oo. We claim that 

(2.13) 0f(r)=O, VA:>2, provided that 0f(r) = r > 1, ;t > 2. 

Note that the growth condition in (12.131 ) is weaker than what is assumed in the 

assumption in Theorem 12. II 

The argument below also applies if is projected on sin^0. First we show that 
= o{r^^) as r ^ oo. Indeed, using < C|x|*^^'+^ we write J^jaije"j(l)j as 

Q(^^-i+T-2-So^ (for some Sq > 0). Let g{t) = (j)-{e'), then from (ITlIl) g{t) satisfies 

g"{t)-eg{t)=h{t), teR 

where 

(2.14) h{t) = 0{e^''-^+^-^>^') t>0. 

Let gi{t) = e^' and ^2(0 = e'" be two fundamental solutions of the homogeneous 
equation, a general solution g{t) is of the form 

g{t)=Clglit)+C2g2it)-^ fg2{s)h{s)ds+^ f g,{s)h{s)ds 

2k Jo 2k Jo 

where ci,C2 are constants. Using (12.141 ) in the above we obtain 

g{t)=c[g,{t)+c',g2{t) + 0{e^'-'+'-'°^'), for? > 1 

where c\,C2 are two constants. Since g{t) = C?(e('^^'+^)') for ? ^ oo, we see that 
Cj = and therefore g{t) = as f — > oo. Equivalently 

(2.15) 0f (r) = 0{r''-^+^-^) i E /. 
With (12.151) we further obtain 

j 
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Consequently (^^{x) = 0(r*^^^+^^^*). Keep doing this for finite steps we obtain 
that 0f decays faster than r^^ at infinity. The asymptotic theory of ODE can be 
similarly used to show that 0f (r) = o(r) as r — >^ 0. 

To get a contradiction, without loss of generality, we may assume that some of 
say 0f (r) > for some r > and 



max , , = max max 



^+ V/iW/ i<j<n\m fj{r) 

By noting ^\{r) = o{r) as r — )• and 0f (r) = o{j) as r — > oo, 0j*//i(r) attains its 
maximum at some point ro € Let w\{r) = <p^{r)/ fi{r). By a direct computa- 
tion, wi(r) satisfies 

(2.16) Aw, +2VW, . ^ + i^w, = i a^je^'^ ( ^^^) . 

Now we apply the maximum principle at r = ro, and obtain 
Awi(ro)<0, and Vwi(ro)=0. 
Since ^ > 1, (12.161 ) yields 



(2.17) £ a,je"^ (ro) < 

because wi (ro) > 0. On the other hand, for j > 2, 

w,(ro)/,(ro) - 0,^(ro) = fMM^ " ^ 

which obviously contradicts ( 12.171 ). Therefore (12.13) is established. When k = I, 
0/ = because by Lemma[2?T] 0/ (r) = CM-(r). By the assumption 0/ (x) = 0(|xp) 
near 0, C = 0. The third statement of Theorem 12. 1 l is established. 

Again by Lemma ITT] the first statement of Theorem 12. H is established. 

Finally, the second statement of Theorem 12.11 is an immediate consequence of 
Lemma 3.1 of [37]. Theorem 12. II is established. □ 



2.2. A Pohozaev identity for global solutions. 
Lemma 2.2. Let u = {ui, ...,Un) be an entire, radial solution of 

-Aui = Y!j=iaije"j , in M?, 



where A is a constant matrix that satisfies (HI). Let 



d = Uj 



m + ^ /^log|T7|(X:a,y^(^))jT7, 

1 /■ " 

a,- = — / e"', = y aijOj, I = min{Zi, ...,/„} 
2n Jr2 
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and 

2k J Br 

Then for some 5q > 

(2.18) =e'''r-''(l+o(r-*)), r>l, 

(2.19) 4£ GiR = £ aijaiRGjR + 2 £ + (9(/?2-'-*) . 

iei ijei iei 

Proof of Lemma |22J 

It is well known that 

(2.20) Ui{x) = I log 1^ - Tj |(£ atje"^^^^)dri + a. 

Indeed, let w, be the function defined by the right hand side of (12.201 ). Then w,- — m,- 
is a hannonic function. Since they both have logarithmic growth at infinity, >v, — 
Ui = c. Evaluating both functions at we have c = Ci. 
Clearly 

(2.21) Ui{x) + li\og\x\ = I log^^f^Vaiy^^^'t/T] +c;. 

271 Jk2 \x\ Y 

Using Zjaije"'^''^ = 0(r-2-*) for some 5o > and r large, we obtain, by elemen- 
tary estimates, 

Ui{r) = -li\ogr + Ci + o{r-^), 

which leads to 

(2.22) Gi = GiR + -^R^''- + OiR"-'---^). 

We arrive at (12.191 ) by using (12.221 ) in the Pohozaev identity for a: 

i ij 

Lemma [Z2l is established. □. 



3. First order estimates 

Let (/jj, be a family of positive, functions on Bi with a uniform bound 

on their positivity and norm: 

(3.1) ^<h'^{x)<C, ||/j?||c3(B,) <C, xeBi, i = \,..,n. 

In the next three sections we consider a sequence of locally defined, fully blown- 
up solutions = (m j , . . . , ) and we shall derive their precise asymptotic behavior 
near their singularity and the precise location of their singularity. Here we abuse the 
notation = {u\,...,u!^^) and it is independent of the one used in the introduction. 
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Specifically we assume that m*^ satisfies the following equation in Bi, the unit 
ball: 

" k 

(3.2) -Au'}=Y,aijh^je"j, i = l,..,n, x€Bi, 

7=1 

with being the only blowup point in Bi : 

maxMf<C(^), V^CCBi\{0}, and 

K Bi 

with bounded oscillation on dBi : 

(3.3) \u'l{x)-u'l{y)\<Cyx,y £dBi, C independent of 
and uniformly bounded energy: 

(3.4) / e"'<C, C is independent of L 

JBi 

Finally we assume that u!^ is a fully blown-up sequence, which means when re- 
scaled according its maximum, {uu} converges to a system of n equations: Let 
u\{0) = maxfij u\ and = and 

(3.5) v\{y) = u'l{e0)-u\{O), y e := B{0,e^'). 

Then v*^ = (vj, ...,v^) converges in Cfg^{M.^) to v = (vi, ..,v„), which satisfies 

f -Avi = 'ljaijhj{0)e''J, R^, i = \,..,n 

(3.6) <^ 

[ /K2e"' <°o, /= 1, vi(0) =0, 

where hj{0) = hmi-^oo/iy(0). 
For the rest of the paper we set 

and m = min{mi, In this section we derive a first order estimate for in 

^k- In l[37il the authors prove that there is a sequence of global radial solutions 
U'' = (C/f , .., of (EH) such that 

(3.7) \u'^{e0)-U^{y)\<C, for b| < ro£,-^ 
From (13.71 ) we have the following spherical Hamack inequahty: 

(3.8) |Mf(£0)-Mf(£^y)| <C 

for all IjI = ly I = < roEj^^ and C is a constant independent of k, r. (13.81) will play 
an essential role in the first order estimate. To improve (13.71 ) we introduce to be 
a harmonic function: 



mi 



(3.9) 
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Obviously 0f (0) = by the mean value theorem and <pj^ is uniformly bounded on 
Si/2 because of (13.31) . Later in section 6, when the results in section 3,4,5 will 
be used to prove the main theorems, the function will be specified when we 
consider the system defined on Riemann surface. 
Let V'' = {V^,..,V^) be the radial solutions of 



(3.10) 



V.^(0) =vf(0), iel 



where vj^ is defined in (13.61 ). It is easy to see that any radial solution V of (13.101 ) 
exists for all r > and e^' € L^(M^). The main result of this section is to prove that 
Vl^{y) + <^f{eky) is the first order approximation to v){y). 

Theorem 3.1. Let A, w*, h!^ = {h\,..,h'^) and be described as above. Then 
for any 5 > 0, there exist ko{5) > 1 and C independent ofk and 5 such that for all 
k > Icq, 

(3.11) |D«(vf(j)-\/.^(y)-0f(£0))| 

( Ce,{l + \y\r"'^'-\-\, m<3, ,-.,-1/2 [«|-oi2 

Definition 3.1. 

\ f k " 

of=— / h\{Q)e^', m\=y aijo], = min{m^, ..,m^}. 
2.K Jr2 

From Theorem \3.1\ it is easy to see that limi;_^.oc = m,. Thus m'^ >2 + S^for some 
5() > independent ofk. 

To prove Theorem 13.11 we have 

(3.12) -^{v\{y) - ^f{e0)) = Y^aijH]{e0)e''^'^y^-^'My\ in a, ( see m) 

1 

where 

(3.13) Hf{-)=h]{-)e^'^\ 

Since 0f (0) = we have //f (0) = /jf (0). Also, the definition of 0f implies that 
— 0^!^(£i-) is a constant ondQ.k. 

To estimate the error term = — 0f (£yf) — Vj^. We find satisfies 

r ^w\{y)+ZjaijH]{e0)e^>w) = -I,-a,y(//}(£0) -//}(0))//, 

(3.14) I 

[ wf(0)=0, /G7, yw\{0) = O{ek), 
where (^f is defined by 

(3.15) e^- = [\*^'-+^'-')^-dt. 
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Since both v'' and V'^ converge to v, = o(l) over any compact subset of M^. 
The first estimate of is the following 

Lemma 3.1. 

(3.16) wf(3;)=o(l)log(l + b|) + (9(l), foryGat. 
Proof: By 

|vf(3')-vf(b|)|<C, VjGH, 
where (r) is the average of on dB, : 

' Inr JdB, ' 

Thus we have e^^-^) = 0(r-2-*) and e^'W = 0(r-2-*) where r = \y\ and 5o > 0. 
Then 

It is easy to use the decay rate of e^'^ , e^'* and the closeness between vf and Vj^ 
to obtain 

r(wf)'(r)=o(l), r>l. 

Hence = o(l)logr and (13.161 ) foUows from this easily. Lemma ITT] is estab- 
lished. □ 

The following estimate is immediately implied by Lemma [TTl 

g4*(.v) < c(l + \y\)-"'+"W for yenk = fi(0,e^"^). 

Before we derive further estimate for we establish a useful estimate for the 
Green's function on Q.^ with respect to the Dirichlet boundary condition: 

Lemma 3.2. Let G{y,ri) be the Green 's function with respect to Dirichlet boundary 
condition on Q.^. For y € Q.^, let 

U = {^^O.k■, lT7|<|3'|/2 } 

E2 = {T7e^^-t; \y-n\<\y\/^ } 

Then in addition for \y\ > 2, 

r C(logb| + |log|T7||), T7Gri, 

(3.17) \G{y,r])-G{0,r])\<l C(logb| + |logb - t]||), t] g £2, 

[ C\y\/\r]\, T]Gr3. 

Proof: The expression for G(j, T]) is 

1 1 Ivl er^y 

G(j,r]) = -^logb-r7| + -log(^|^-77|), y,^ ^Q.^. 

In particular 

G(0,77) = -;^log|77l + :^logei"\ t? e ^^-t- 
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Therefore we write G{y,ri) — 0(0,7]) as 

1 



(3.18) G(3;,77)-G(0,77) 



-log^ + J-logl^-^ 



The proof of (|3.17| i for T] G Zj is obvious. For Tj G ( 13.171) also obviously 
holds if either IjI or |t] | is less than because in this case 



log 



bl 



■£l\y\n\ 



<c. 



Consequently 



|G(j,T])-G(0,T7)| < c(log|T]| + |logb-T7||+C) 
< C(log|T]| + |log|3;-T7||). 

Therefore for 77 € we only need to consider the case when \y\,\^ 
this case it is immediate to observe that 

y 



> is, 



7^-1 



In 



logl^-e^^blr?! 



where Z(-, •) is the angle between two unit vectors. Thus for T] G £2 we only con- 
sider the situation when b|,|T]| > |£^\ Z(||y,^) < |. For this case we estimate 
G(j, 77 )-G(0, 77) as follows: 

|G(3^,T])-G(0,T])|<|G(3;,T])| + |G(0,T])| 

|G(0,T])| <Clog|3;|. 



\GiyM < 



1 I, t ,1 1 , 8 1 
-llogb-r]|| + ^log- + - 



271 



1 



riW 



< C(log|3;| + |log|3;-77||) 
where the last inequality holds because 



\y-ri\ < 



' bP 



■'7l<cbl, 



which implies 



h y 



■T]|| <C(|log|3;-T]||+log| 



The second case of (13.171 ) (when rj G £2) is proved. 

For T] G £3, we first consider when |t] | > 2[3'|. In this case 

bl 



|log- 



|l0£ 



\ri-y\ 

For the second term, since ri,y £ D.k and |t] | > 213^1, we have bl 1^ I < \^k'^' conse- 
quently 



bh 



bl 
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So (13.171 ) is proved in this case. Now we consider '-^ < |t7 | < 2|y| and ji] — j| > -y . 



For the first term we have 



log- 



\y-ri 

For the second term, we want to show 



\y\ 



(3.19) 



log 



bl 



,-2 



<c<c 



\y[ 
1^1 



If either \y\ < jle^ ^ or |t] | < jle^ ^ we have 



,y_ 
M 



£,-2' - 16' 



therefore (13.191 ) obviously holds. For ||e^ ^ < \y\,\ri\ <e^ \ using Ij — t]| > \ \y\ 
we obtain easily 



y_ b|T] 

bl 



-2 



> 



Therefore (13.171 ) is proved in all cases. Lemma |X2] is established. □ 

Proof of Theorem I3.lt First we prove (13.111 ) for a = 0. We consider the case 
m < 3, the proof for the case m > 3 is similar. By way of contradiction, we assume 

max/e/|wf(j)| 

Ak := max — - — , , J , , > °°- 

ye^i,Ek{\ + \y\f+^-"' 

Suppose A^: is attained at yk G O-t for some /q G /. We thus define 

_k, . _ w\{y) 

Here we require 5 to be small so that m — 2 — d > (Thus 3 — m + 5 < 1). It 
follows from the definition of Aj^ that for y £ Q.k 



(3.20) \w'^iy)\ 



\w-{y) 



\3+S- 



< 



:i+bi) 



3+5- 



A,4(i + blP+'-'"(i + b/tl)'+^-'" - (1 + ^1)'+' 

(1 + 



The equation for is 

(3.21) -A^'^(y) = Y^aijh'j{0)e^'w) + o{l)- 



\1— m 



(1 + b.l) 



3+5— m ' 



for i G /. Here is given by (13.151) . ^j^ converges to v; in Cfg^{M.-^). Besides, 
we also have (0) = for all / and Vw\{0) = o(l). If a subsequence of yt stays 
bounded, then along a subsequence vv*^ = (vvj, converges to vv = (wi, ...,vv„) 
that satisfies 



-Awi = l^j aijhj {0)e"j wj, , / G /, 

w,-(o)=o, Vwi(o)=o, \wi{y)\<c{i + \y\y+^-'\ 



y^- 
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Thanks to (1) of Theorem 12. II 

Since Vvvi (0) = we have c\=C2 = 0, thus w; = for all /. On the other hand, the 
fact that w\jjk) = ±1 for some /q £ ^ implies that = ±1 where y is the limit 

of j^. This contradiction means that oo. Next we shall show a contradiction if 
\yk\ °°- By the Green's representation formula for vv^, 

wKy)= I G{y,^){-^w'l{^))d^+w\\^a, 

where vvf l^f^* is the boundary value of vv, on dO-u ( which is a constant). From 
(IHOl ) and (feB we have 



(i + b/tl)^+^-"' A;t(i + bi|)3+^-'"' 

Thus for some i ^ I we have 
(3.22) \ = \wHyk)-wUO)\ 

< C j^^ |G(,„r7) - G(0,r7)| ^ ^ |^^|^3^,_„, + ^^^^ ^ b.|)3+^- J ' 

where the constant on the boundary is canceled out. To compute the right hand 
side of the above, we decompose the Q.k as = Zi U Z2 U £3 as in Lemma 13.21 
Using (13.171 ) we have 

/ |G(j,,r7)-G(0,r7)|(l + |r7|)3+^-2'"jr] = O(l)(logb,|)(l + b,|)(^+^-2'«)^ 
where 

" (1 + b-tl)", «>o, 

(i+bir" = <( iog(i+b^i), «=o, 

1, a<0. 



\G{yu,n) - G(0, T7)|(l + |T7 1)3+5-2.^,7 = o{l){\ + b^D^+^-^M. 
Hence 



2— m 



1 i^f"'") - "'"■■'>i atMP-.> '' = OCX' + 

Similarly we can compute the other term: 

y^jG(,,,r7)-G(0,r7)|-l^±M^ 

By the computations above we see that the right hand side of ( 13.221 ) is o(l), a 
contradiction to the left hand side of (I3.221 i. Thus (13.1 11 1 is established for a = 
0. The estimates for |a| = 1 and 2 follow easily by scaling and standard elliptic 
estimates. Therefore Theorem 13. H is completely proved. □ 
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4. Second order estimates 

In this section we improve the estimates in Theorem 13. II for m < 4 and m = 4, 
respectively. Let pij^ be the maximum point of vf(-) — 0f (£<;■). The following 
lemma estimates the location of pi^k- 

Lemma 4.1. pijt^ = 0{ek), i € /. 

Proof : Applying Theorem ITT] to vf — 0f (5-) on Bi: 

(4.1) Z)«(vf(3;)-0f(£0)) 

= D^{vt{\y\)) + 0{e,), \y\<\, |a|= 0,1,2. 

The equation for V^^ is 

iVh"ir) + -iVtnr) + t aijH]{0)e''nr) = 0, r > 0. 

From {V^y{Q) = 0, we see hm,^o(^-'')'(r)/r = (V,.*^)"(0). Thus 

(4.2) {vt)"{Q) = -l-t «/;^}(0).^' < -C 

for some C > independent of k. Since pi^u is the maximum point of v^(-) — 
(£yf)' we deduce from (1411 ) that (V^*')'(|;?,,yt|) = 0(ei), thus from (l42l ) we have 
/'i./t = 0{ek). Lemma l4~T] is established. □ 

The main result in this section is to find the eu approximation to vf (•) — 0f (e^;-). 
It is most convenient to write the expansion around one of the pi,k^. We choose 
p\^k and shall use <I>* = (<I>j, ..,<I>^) to denote the projection of v^(-) — 0f (£*:•) onto 
i'/7a«{sin 0,cos 0}. i.e. 

(4.3) <I>f(rcos0,rsin0) =ei(G^;(r) cos + G^;(r) sine), /G/ 

with Gf j(r) {t = 1,2) satisfying some ordinary differential equations to be specified 
later. 

Set v^'* = (v|'*^, ..,v^'*) as 

(4.4) v]H-) = v\{.+pu)-<i>h£k{- + Pi,k)) 
in 

(4.5) D.l,k■■={^■,^+p\,k&^k}■ 

Using Vvf (0) = 0{ek) (by TheoremllD and 0f (0) = we observe that 

(4.6) v]^\0)=v'l{pu)-^h£kPu) 

= vf (0) + Vvf (0) ■ pi,k + + 0{pi,k£k) 

= v\{0) + O{el) 
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The equation that v' *^ satisfies is (combining (13.111 ) and (13.121) ) 



(4.7) 



^v]'' + 'LUaijHf'{y)e'/ =0, in Hi,, 



[Vv;''(0)=0, Vv]^'{0) = O{e,), i = 2,...,n 

where H^'' = {hI'\ is defined by (see (I3l3]) ) 

(4.8) hI'\-) = H^isk ■ +ekpu) = h^Sk ■ +ekPu)e 
Trivially 



(4.9) 



Hl\0)=hm + O{et). 



In the coordinate around p\ k and we seek to approximate v]''^ in flu. The first 
term in the approximation of v^ '^ is V'^. Here we note that the domain ^ is shifted 
from the ball Hj. by p\ k- 

We shall use five steps to establish an approximation of v]'^ without distinguish- 
ing m = 4 or not. 



Step one: 



Let w''*^ = (wj *^, .., w^'*^) be the difference between v*^ and V 



k. 



^]\y) = Vi\y)-yn\y\). y&^it- 

Taking the difference between (I4.7l i and (13.101 ). we have 

j j 

which is 

Aw^^+I«,/.5(0)//(^^^-l) = a 

Here we observe that the oscillation of on dQ.)r is O(e^). Indeed, recall that Hi ,t 
is the shift of the large ball Q.^ by pi^. Let y\,y2 & dQ.i^k, one can find yj, such that 

bsl = \y2\ and 1^3 —yi\< Ce^. Since (V^-*^)'(f) ~ r^^ for r > 1 and \yi \ ~ we 
have 

(4.10) vi'iyi)-vHy2) = vHyi)-v,Hy3) = o{el). 

With (14.101 ) we further write the equation for w- as 



\.k 



(4.11) 



w'f\Q) = 0{el), i G /, w\'' = 0{£l) on da.^u, 



[ Vw;'*^(0)=0, Ww\^\Q)=0{et), /G/\{1} 



\,k 
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where 



Similar to Theorem 13. ll we also have 

Lemma 4.2. For any 5 > 0, there exists kQ{5) > 1 such that for some C > inde- 
pendent ofk and 5, the following estimate holds for all k > k^: 

Proof: Using the definition of v' *^ and Theorem 13. ll we have 

On the other hand we clearly have 

\vhy)-Vhy + Pu)\<Ceu{\ + \y\)-' 
by mean value theorem and the estimate of VVj^. Lemma l42l is established. □ 
Using Lemma |42l and ( I4.9I ) we now rewrite , clearly 

£f = -l^ai/'{H]''{y) - h]{0)) + {H^{y) - h^m^]'" + 0{{w]''f). 
j 

By Lemma |42l and ( I4.9I ). the last two terms are 0{e^{\ + l)'])^^'") regardless 
whether m > 3 or not. Thus 

(4.14) Ef = -l^a^je^f{H'/{y)-h){0)) + O{e^){l + \y\f-'" 
j 

= -l^a,je''f{Hf\y)-H^\0)) + O{el){l + \y\f-'" 
j 

where in the last step we used (14.91 ) again. 
Step Two: Estimate of the radial part of *^: 

k,0 

h 



Let g*^'" = {g{ , ...,gn ) be the radial part of 



1 r"^^ 

si {^) = IT' I ^i' {rco&d,rsind)dd. 
2n Jo 



(4.15) 



Due to the radial symmetry of v!^, g^'^ satisfies 

[gfiO) = Oiei), iel, ig\\0)=0. 
where ( for simplicity we omit k in L,) 

T k.Q k,0 , 1 d 1,0 r-i ,/t/n\ V* kfl 
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We claim that for 5 > 0, there exists ^o(5) > 1 such that for all k > kg 

(4.16) \g':-\r)\<CeUl+rf-'"^', 0<r<e,' 

holds for some C independent of k and d. So, gf'*'(r) can be discarded as an error 
term. 

To prove ( 14.161 ). we first observe that m < 4 and by ( 14.151 ) 

<Ce2(l+r)2-'"+5/2. 

Let = {f\^..,fi^) be the solution of 



dr 



fk \ \_d_fk T k,0 

■Ji ' r drJi ~ ^'6 ' 



Then elementary estimate shows 

(4.17) \ff{r)\<Cel{\+rf-'"+'. 
If m > |, we claim 

(4.18) \gfi^r)-ff{r)\<Cel(\+rf, 0<r<e,\ 
Indeed, let = g^-^ — f^, then clearly 



k 




dr 



ff(0)=0, 



where 



j 



By considering ^ + as A in and g'^ a solution with boundary oscillation 
in B(0,e^'), we obtain ( 14.181 ) by the argument for Theorem 13.11 Here we note 
that to apply Theorem 12.11 it is essential to require 6 — 2m + 5 < 1 (it holds if 
m > 5/2), the estimate on the Green's function in Lemma 1X2] and the condition 
^g^(O) = 0. Since m < 4, 0{el){\ + rf is part of the error. 

If m < I we apply the same ideas by adding more correction functions to g*^ ": 
Let A'^ satisfy 2 + (2 — m)N < 1, we add N correcting functions to make the right 
hand side of the equation of the order C?(£^(l + j-)(2-»i)^+5^ Note that each cor- 
rection can be discarded as an error in the sense that they are smaller than the right 
hand side of (14.471) . Using 2 + 2N — mN < 1 and the argument in the proof of 
Theorem 13. II we obtain (14.161 ). 

Step Three: The projection on sin 6 and cos 6 

In this step we consider the projection of w^ *^ over cos 6 and sin 0, respectively: 

£kG'u{r) = — w \r,d)cosddd, £,G^,.(r) = — / w ''(r, 0) sin 

271 Jo ' 271 Jo 

Let 

(4. 19) of = ekG\i{r) cos G + e^G^ ,(r) sin G 
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clearly G\ , and G2 , solve the following linear systems: For < r < and t = 1,2 

,d^^ld 1 , 
' dr^ r dr r^' 



(4-20) ^^S + Laim^'^, 



J 



= -Y^a^jd,H]-\0)re''^ +o{ek){l+r 
Then O'' solves 



2— m 

J 

Uijllj\\J jt^ ' 

J 



(4.21) A<I>f + £a,y/j5(0)/.<I>5 

= -et'£a^j{d,Hf\0)y:+d2H]\0)y2)e''f + O{el){l + \ 



\2— m 



By the long behavior of w' *^ (Lemma 14^21 ) we have 

(4.22, iG!,,Mi+iGLMi<(^i:!:!r'''"-'' 



C(l+r)^ m>3. 
Note that Gf ,(0) = and Gf ,.(r) = 0(r) near 0. 

Step Four: Projection of onto higher frequencies 

Let 

be the projection of w^ *^ on sin/6. In this step we first establish a preliminary 
estimate for all these projections: 

Lemma 4.3. There exist lo>3 and C > independent ofk, I such that 
\gY{r)\<Celr\ 0<r<e,\ W > /q- 

Proof: 

By (14. 111 ). (14.141 ) and Lemma l4!2lg^'^ satisfies 



(4.23) 



I g-''(o) = o. 



where, applying the Taylor expansion of H'-^ up to the second order, 

(4.24) \hiki{r)\<Cel{\+rf-'"+^l^ 
for some C independent of k and /. Thus 

(4.25) Agf'-^gf' + £a,,-/^5(0).^gf > -c,el{l + rf-^' 

for some co > 0. Let 

g{r) = - t^{l+sf-"'+^ds + co— s\l+sf-"'+^ds + r\ 
4 Jr s 4 Jo 
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Then clearly g{r) > for r > 0, g solves 

f / + ^g'-4g(r) = -cor2(l+r)-'"+^ r > 0, 

(4.26) <^ 

I g(0)=g'(0) = 0. 

and 

g(r) = r2 + 0(/-'«+^), r>l, g{r) <Crhog{- + \), r<l. 

r 

Clearly by Lemma 14^2] 

(4.27) elg{e,') > max |wj'^(y)| > {g'/ {e^')\. 

The reason that we include in the definition of g{-) is because by Lemma l4T2] we 
only know w'r\x) = 0(£,(1 + |x|)^) for m > 3. Let g' = {g\,..,g',) = fe, ..,g), 
then it is easy to see that for Iq sufficiently large and / > Iq 

(4.28) Agf-^gf + £a,;/j5(0)^'^g5 

= -elcoil + rr-^' + j:a,h)iO)e'}g) - '^g^ 

< -£fco(l + r)2-"'+^ 
To prove Lemma 1431 it is enough to show 

(4.29) \g-''{r)\<Cgl 0<r<e^\ / > /q. 

with C independent of k and /. To this end, we shall use ( 14.251 ) and the initial value 

(4.30) =o(r) near 0. 

In the following we use the argument in the proof of Theorem 12.11 Let \{/f = 
^i-gY, by ^^,^M,^M and (l430l) 

' ^Wf-'^wf + L; aijh) (0)// v^} < 0, < r < \ 

< 

(r) = o(r) near 0, Yii^k^) > ^ 

Our goal is to prove \l/^{r) > for all < r < To do this, let f'' = (/f , ...,/^) 
be the positive solution to the homogeneous system: 

such that 

/;^(r) ~ r, near 0, /;^(r) ~ - near oo, /;?'(r) > for all r > 0. 
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If ^ff is not always non-negative we assume 

. W'tir) . . wHr) 
mm , = mm mm -LL^—L < 0. 

R+ /[(r) !<!<« R+ f!'{r) 

Suppose the minimum is attained at r^. By the behavior of i//f and f^, we have 
< < £^"^ Let w\ = V^f//f , then w\ satisfies 

Evaluating both sides at r^, the left hand side is strictly positive while the right 
hand side is non-positive by the definition of w\. This contradiction proves (|4.29l l. 

Finally since g{r) = 0{r^logj) near 0, (r) = 0{e^r^) near for / > Zq > 3. 
Lemma |43] is established. □ 

Lemma 4.4. Given 5 > 0, there exist C{5) > independent ofk,l and ko{5) > 1 
such that for / > 3 and k>kQ 

(4.31) \g^^'(r)\<Cet-'-'^\e,rU^ey{l+rf-'"+^ r <\e,' 
where 

m — 2 — 5 if m<3 



I o s\+ \ m-l-b if m'_ 
{m-2-by=< ^ ' 

y \ — o ij m>i. 



Proof: 

By Lemma 1431 (14.231 ) can be rewritten as as 

^d?^-rd-r-7^^'' =^'-^'' (^)=° 

where 

(4.32) \kM{r)\<Cei{l+rf-'"+^, k>ko{5) 
for some C independent of k and /.By standard ODE theory 

(4.33) gf'{r)=c[„r' + c'2„r-'-^^ p s-'+%,,{s)ds - ^— s'+'hki{s)ds. 

When r ^ 0, it is easy to see that the last two terms in (14.331 ) both tend to 0. Thus 
'^2ki ~ ^- ^3 ^4 represent the last two terms in (|4.33l l. respectively. By (14.321 ). 
we have 



(4.34) |/3(,)| + |/4(,)|<^(i+,)4-+5 



where C is independent of k and Z. On the other hand, using the information of w 
for r ~ we know 

(4.35) \gf'(r)\<Cet-'-'^\ r^e,\ 
In regard to (14.341 ) and (14.331 ) we have 

(4.36) \c\u\<Cet^''^'^'- 
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Hence (14.311) follows immediately. Lemma l4!4l is established. □ 

Let z*^ = {z\, ...jZ^) be the projection of w^'* to span{smW,cosW,l > 2}, i.e. 

4 = l(^-''Wsin/e+f^'(r)cos/0' 



1=2 



where g^/ is the projection of w-'*^ on cos/0, g'^'' has similar estimates as that for 
g^/. Then Lemma |44l leads to 



(4.37) \zHy)\ < C£f (1 + bl)'-'"+^ bl < 
for m <3. However for m> 3, Lemma |44] only gives 

(4.38) \,^(y)\<Cel-'{l + \y\f+Celil + \y\r-'"+', \y\ < 



In the following we shall get rid of the first term on the right hand side of (14.381 ). 
To this end, we need to evaluate the value of w^ *^ on d^i^k- 

Lemma 4.5. (a) Ifm<4, then 

(4.39) \w[\y)\ <Cef -Hoge^\ y e da^^f 
(b) Ifm = 4 and \m\ — A\ < C /\oge^^ for all i, then 

(4.40) \w]'\y)\<Cel{\ogekf, J G 5^1,^ 

Remark 4.1. The assumption — 4| < C/log e^^ when m = 4is natural and will 
be justified in the proof of the main theorems in section 6. We also remark that in 
i\4.391l we use £™ ^^log 1/e^ instead of the crude £^^^^^ as before. 

Proof: 

The Green's representation formula for w]''^ gives (see (14.111 )) 
O{el)=w]\0)=j G{0,n){'£a,jh]\0)e'''w]\r^)-E^{Tl))dn 

(4.41) - / ^yG{0,^)w'r\^)dS^ 

where G{-, •) is the Green's function on with respect to the Dirichlet boundary 
condition. We observe that 

(4.42) G(0, T] ) = - log I T] I + log 1 + 0{£i) . 

Indeed, since D-i ^ is a translation of 11^ by pi j^ ( recall pi j^ = 0{£k)), the oscillation 
of log [t] I on dQ.\ u is 0{e^). Thus the oscillation of the regular part of G(0, T]) is 
O^el), which leads to (14.421 ). On the other hand let w- be the average of w- on 

^Hijt, using the fact that v^- is constant on ^^i^^ and Vj^ has oscillation O(e^) on 
dD.i kwe have 

(4.43) w]-\y)=w]'' + Oi£l), G ^^i,,. 
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-W; 



Thus 

(4.44) -/ dyG{0,ri)w'r\ri)dS,^=wj^' + O{ei). 

By using (14.421) . (14.431 ) and (14.441 ) in (14.411 ) we have 

+ = |^^^(i-log^)(£«,,./^5(0)//vv;.'^(r7)-£f(r7))c/r7 

To evaluate the right hand side, we divide Hijt into a symmetric part: D[ := 
B(0,e^' — |/?i,/t|) and a nonsymmetric part: Hi ,t \Di and use /i and I2 to rep- 
resent the corresponding integrals on them. If T] G Hi ;t it is easy to see from 
\pi,k\ = 0{ek) that 

£-1 

log(^) = 0{e^) and iH^A^il = 0(1). 
Moreover for 77 € Hi \Di, by (14.131 ) and (14.141 ) 

w;Ar]) = 0(£r^), £f(r]) = 0(£r*). 

2\ 



Combining these facts we have I2 = 0{e^ 
To evaluate /j, since log(£^^'/|T] |) is 
of wY and onto 1 remain. By (14.161 ) and (14.141 ) 



To evaluate /j, since log(e^ Vl^l) radial, by symmetry only the projections 

i,k 



J D[ 



(Note that we use the fact that f'^^W < C(l + jjl)"™* for C independent of L) There- 
fore if m < 4, /i = 0{ef log e,^^. If m = 4 and |mf -4| < Cloge^"' we have 
(1 + r)^'" <C{\+r) , thus by elementary computation 

\h\<Cel{\ogetf. 

Lemma l43] is established. □ 

By Lemma 1431 (14.351 ) can be replaced by 

\gY{r)\<Ce'r^-\ r~£-i. 
Correspondingly, the estimate for Cj^, becomes 

\C\kl \ ^ 

Then it is easy to see that the first term in (14.381 ) can be removed. 
Step five: 

For the projection on s/7<3«{sin 0,cos 0}, we write (14.201 ) as 

(^+~ ->!,■=%). o<.<s- 

where h is the collection of other terms. By (14.221) 

\h{r)\<C{\+rY-'"+^. 
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Then 



(4.45) G\ i = cur + — - x / h{s)ds - — / s^h{s)ds. 

r 2 Jr 2 Jo 

Since Gij is bounded near 0, cjk = 0. Using G\ ,(e^') = 0{e^^^^^) we have 

I _ I ^ ^„m-2-S 

\c\k\ < Ce^. 
Then it is easy to see from (14.451 ) that 

(4.46) \Gli{r)\ < Cr{l + r)^-'"+\ t = 1,2, whenm > 3. 

Similarly (14.461) also holds for G\ 

Combining the results in the five steps we arrive at the following estimate with- 
out distinguishing m < 4 or not. 



Theorem 4.1. Given 5 > 0, there exist C{5) >0, kQ{5) > 1 such that for \y\ < 
£^7 ^nd \a\ = 0, 1, the following holds for all k>ko 

(4.47) \D"{vl\y)-VHy)-^Hy))\ < Cei{l + bD^—l"^^ 
where 

of (3^) = ek{G\jir) cos d + G^ir) sin 0) 

with 

(4.48) I^mWI < Cr(l + r)2-'«+^ t = 1,2. 

Note that the estimate for [a| =0 follows directly from the five steps. The 
estimate for | a | = 1 follows from standard gradient estimate for elliptic equations. 

Theorem 14.11 does not distinguish m < 4 or m = 4. In the following we apply 
Theorem 14. ll to obtain more accurate estimates for m < 4 and m = 4, respectively. 
Both results in the sequel (Theorem 14.21 and Theorem I4.3l l play a crucial role in 
determining the location of maximum points of bubbling solutions. 

Theorem 4.2. Suppose m < 4, then for \y\ <e^^/2 and i € / 

(4.49) |D« (vf (y) - 0f (£0) - V^iy - p,^u) -^{y- Pii)) I 
< C£2(l + |3;|)4-'«*-'log(2+|3;|) |a|=0,l 

where vf , 0f , V^*',<I>f are defined by <\3.5\) . ( I J. 91 ). liS.lOll and ( 14.51 ). respectively. More- 
over Gf ; (t = 1,2,/ G /) satisfy 

(4.50) |Gf,.(r)| <Cr(l+r)2-™' 0<r<£^-i. 
Proof: 
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We use the same notations as in the proof of Proposition 14.11 First we consider 
the radial part of w'*^: Recall m < 4. Using Proposition 14. 1 l and (14.151 ) we have 

= -£a,./.5(0)//gf-|£a,,A//|^0)r2// + O(£,2)(l + r)-'" 

j j 

= 0{el){\+rf-'"\ 0<r<£^-i, 

and 

Multiplying r on both sides of the equation and integrating, we obtain 

\gf\r)\<Cei\og{2 + r). 
Next we consider the projection of w^'* on sin 6 and cos 6. Let G)*^ = (g)^ G),^) 



be 

Then co'' satisfies 



k l,k k,0 

(Oi=wr -g. . 



Acof + £ aijhp)e''jco'; = E^, Hi,- 

where E^^ is the projection of the original right hand side on the subspace spanned 

by sin^0 and cosku {k = 1,2, ..). Since gj^'" is a radial function, from the asymptotic 
behavior of g*^ " it is easy to see that the oscillation of it on dQ.i^k is 0{e^), therefore 
the oscillation of cof on dO-i^i^ is 0{e^). Let cbf be the average of (of on dQ.i^k, then 
by G)f (0) = we have 

Using (14.421 ) in the equation above we have = 0{el), thus 

(4.51) «f = 0(£f), on 5^1,,. 

Since we have known that £i;Gf,(£^') = 0(£™*^^), we can improve the estimate 
of using (14.451) . The estimate of Gf now is 

\GUr)\<Cr{\+rf-'"\ r < e^' 

and 

\j-Gli\<Cr{l+r)'-'"\ 0<r<e^'/2 

which leads to 

(4.52) |V<I>f(j)| < Cetr{\+r)'-'"\ \y\ < e,'/2. 
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As far as the projection of w^ '^ on higher frequencies is concerned, since we have 
(14.5 lb . for g'^'' we now have, instead of Lemma l4.4l 

(4.53) \gY{r)\ < Cef-\e,r)' + ^^e^ [l + rf-'-\l > 2. 

As before we let z| = wj''^ — <I>f , then z'} satisfies 

Az'^ = Oiei)\ym\ + \y\)-'"\ auCR\ 

= 0(£2) on d^i^k- 
Because of (14.531 ) we have 

\zHy)\<Cel{l + \y\)'-'"\ \y\<e,\ 
By standard re-scahng method 

(4.54) \VzHy)\ < Ceiil + \y\)'-"'\ \y\ < e^' /2. 
We have established 

\D"{v]^\y)-Vny)-M{y)\ < Cel{\ + b|)^-'"*-l«llog(2+ b|) 

for |a| = 0, 1 and \y\ < e^'/2. Recall that v- '*^ is defined in (14.41 ). Instead of using 
the coordinate around pi^k we use the coordinate around the origin to obtain (14.49) . 
Theorem I4.2l is established. □ 



Theorem 4.3. Ifm = 4 and \m\ —A\ < C/log£^ ' for all i € /, then we have, for 
\y\ <£;t V2 and i el 

(4.55) \D'=^^v\{y)-<i,^{euy)-VHy-Pu)-^'{y-Pu))\ 
< Ce2(i + |3;|)-l«l(iog(2+bl))2. |a|=0,l, 

where is of the form stated in ( 14.51 ) with Gf = l,2j satisfying 

(4.56) |Gf,.(r)| <Cr(l + r)-2, < r < £^-\ iel. 
Proof: 

For m = 4 we use (14.161) to write (14.151 ) as 

< 

Note that (1 <C{\ + r) . Multiplying both sides by r and integrating, we 

have 

8i\r) = -e^l^au^^[t'e'''^'\\ogr-\ogt)dt + 0{e^,)\og{l+r). 

To evaluate the integral, we use (I2.21l i to write 

///'*^(0)e^*W = e'h-'"' + Oit-^-"^), t > 1 
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for some £o > 0- Under the assumption of Theorem 14. 31 \m!^ —A\ < C/log£^' we 
have t-'"' = 0{t-^), thus 

\g]'\r)\<Cel{\og{2 + r))\ < r < e,\ 

The projection on higher frequencies has the same estimates as in the case for 
m < 4. Specifically, let cof^ be the same as in Theorem 14.21 Then (14.511 ) also holds. 
Correspondingly (|4.52| i. (14.531 ) and (|4.54| i still hold with m*^ = 4 by the same proof. 
Theorem l4.3l is established. □ 

5. Location of the blowup points 

In this section we determine the locations of the blowup points in Theorem 
and Theorem 



Theorem 5.1. Under the same assumptions as in Theorem 

(5-1) \L[^§^ + MhO))of\<Cef-\ / = 1,2, 

where C is independent ofk. On the other hand, under the assumptions in Theorem 
]we have 



(5-2) \L{^§^ + di<t>hO)yh<Ceiloge^\ 1 = 1,2, 

where of is defined in Definition \3.1\ 

Proof of Theorem ISHl Recall Hf{Eky) = h'^{£ky)e'>'-^'"^y^ and ///'^ is defined in 
KE\i . Let Q.k = B {0,8,^^/2), we use the following Pohozaev identity for the equa- 
tion for v^'*^: For <^ eS\ 

(5.3) £ / d^Hl-'(y}e^"^^ 

According to the definition of ///'^ in (14.81) 

d^Hl\y) = eud^HfiO) + 1 eld^iHf{0)y' + 0{el{\ + \y\f). 

i=\ 

Using the expansion of v- '^ in Proposition 14. 11 ( which holds for m < 4 and m = 4) 
we have 



{eud^Hf{Q) + f eld^,Hf{0)y^ + 0{el{\ + \y\f)) 
1=1 

.(^e'''^y\l+^'}) + Oiei)il + \y\)^-^"'+'yy. 
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By symmetry we have 

2 

Also by elementary estimates we have 

/ {teld^iH^{0)y' + O{ei{l + \y\f)){e''-(y^^^ + O{el){l + \^^^^^^ 
and 

/ e,d^Ht{0)O{ei){l + \y\)'-^'"+' = 0{ef-'). 

Thus 

= [ £,5^//f(0)^>^*W + O(£|(l + b|)2)e^*W. 

The only difference on whether m < 4 or m = 4 is on this term. For m < 4, direct 
computation gives 

/ d^Hl\y)e"^y^ 
= SkdtHfiQ) [ e'^-^yUy + Oief-'). 

On the other hand for m = 4, using the closeness between and 4, we have 

/ d^Hl'\y)e^''''^y^ 

= iKEkd^ (log //f ) (0) of + 0{el log 1 ) . 

Next we estimate the right hand side of the Pohozaev identity. The key point in 
this part is that the only difference between Theorem 14. 2 l and Theorem l4.3l is on the 
radial part of w-'^: gl'^ . The projections of w- on higher frequencies {cof ) have 
the same estimates, i.e. (14.511 )- (14.541 ) hold for m < 4 and m = 4. The detail is as 
follows: First by the decay rate of v\'^ one sees easily that 

/ {^■v)l^e^^'Hl-' = 0{ef-'). 

Before evaluating the remaining two terms we first observe that by (14.521) and 
(14.541) that 
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Moreover, by symmetry we have 
Thus 

= E/- o"((av«'+sf) + o(e;"'-'))p5(v?+sf) + o(e;"'-')) 

= 0{ef-'). 
Similarly in the last integral: 

i j ■JdQ.t 

Theorem 15 -H is established. □ 



6. The leading term for p'^ ^ p ^Ti 

In this section we complete the proofs of Theorem II .21 Theorem 1 1.31 and Theo- 
rem [TTT] The notation refer to the one in the introduction. Let 

(6.1) @'l = u\- log / hie"-dVg, i G /. 

Jm 

Then we have 

(6.2) - Aj0f = £ aijp](hje'^' - 1). 

7=1 

where 

(6.3) //i,Ay„ = l, Vol{M) = l. 

JM 

is a sequence of blow-up solutions. We use pk to denote the point 
where the maximum of 0*^ on M is taken. Take the local coordinates around pu- 
Then ds^ has the form e^^^n)i^dy\ + dy2) where 

I Vva(O) I = 0, v^(0) =0, ^^lf= -2Ke^, 

K is the Gauss curvature. In local coordinates, (|6.2| i becomes 

(6.4) -A0f=£fl,-,p}e'^(V®*-l), in Bg, 
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for some 5o small. Let fj^ solve 

' AfHLUaijp'^eV, Bs, 

^ f^{o) = o, v/;^(o) = o. 



(6.5) 

and fj^ has bounded oscillation on dBg^. Set 



0f = 0;(w + -)-/;' and ^ = p^hi{pk + ■)e'^^^'' 
then satisfies 

(6.6) 



Note that from the definition of h^, ff and ^f, the following can be verified by direct 
computation: 



Ah'^iO)=pf(Ahi{pk) + hi{pk)i-2K{pk) + ZU''iJPj))- 



(6.7) 



Let Mk = max,£/ maxM , = e 2^'', 
(6.8) vf (y) = (e,y) +21oge, = 0f(e,y) -M,, 

then we have 



(6.9) 



Av1(y) + i: atjh%eky)e'^J^^ = 0, b| < ^£,"1. 



The following lemma proves that the ©* is a fully blown up sequence. 

Lemma 6.1. Along a subsequence (vj, ...,t^) converge in Cf^^{M?) to a system of 
n equations. 

Proof: By way of contradiction we assume that only I {I < n) components of 
(vj,..,v^) converge to a system of / equations. Without loss of generahty we as- 
sume that the first / components of ( v j , . . , ) converge to ( v 1 , . . , v/ ) that satisfies 



-Av,- = ^ aijhje^' , ? = 1 ,..,/. i 



in J 



7=1 



where hi = Pihi{p). Let 

Then the entire solution (vi, ...,v;) with finite energy satisfies 
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Let J = {1, ...,Z}, then by Theorem C in ll38l (vi, ...,v;) also satisfies 
1 



By (163 



2;r 271 JB(pk.,S) Jbs 
Thus by letting ^ oo we have 

^ > Cm-, /=!,..,/. 

Let = ^ — Ci.v, ' G J, then easy to see 

(6.10) -L(A^(p)_A^(a.,)) = -2^(£a,ya;>-2>, -£a;y^;^y <0. 

ieJ jeJ Lj 

Since A7(av) = 0, (16.101 ) is a violation of the definition of Fi in the introduction. 
Lemma l6?T] is established. □ 

Let (j)^ = be the harmonic function that takes at and that kills 

the oscillation of ©'^ on dBg^- The first term in the expansion of — 0^(£yf) is 
U^{--P\,k) where U'' = (t/f , satisfies 

\ t/f(0) = vt(0), /G/. 
Moreover, + \o%{p'[hi{pk)) satisfies 

-A([/f + log(pf = £a,■,■.^^'°^(''^^(^*» Ml 

Set 

(6.12) 0^ = ^ [ P^hi{p,)e^', m\ = f^a,ja]. 
Correspondingly 

= min{mi , m=limm^. 

Later we shall show that m is equal to the one defined in (11.61) . In regard to Lemma 
I2.2l we set as 

(6. 13) c'l = f/f (0) + log(pf/i,(w)) + ^ log 1^ 1 1 aijh%0)e'''d7i 

If" * 
= ©f (pO -M^ + log(pf/j;(;7^)) + — y^^log |t]| £ aijp';hj{pt)e^jdTl, 



and 

(6.14) c, = limcf. 
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According to Proposition 14.11 the expansion of can be written as 

(6.15) vf (3;) = Uj^iy - p,,u) + ^h^kj) + ^{y - Pu) + 0{el){\ + b|)^-'"+'= 
on 6(0,^)6^'), where pyj^ = 0{ek), 

^'^{y) = £,(Gf,.(r) cos + Gl,(r) sine) 

such that 

\Gl,{r)\<Cr{l+rf-'"\ 0<r<doe^\ f = 1,2. 
Here we note that the Green's function is of the form 

) = - ^ log |x - T] + y{x, T] ) 

where ;^ is a cut-off function such that ;^ = 1 in B^,^ for some £1 > {2e\ is less 
than the injectivity radius of M) and % = outside B2ei • the sequel we always 
assume So < £i/4. 

The Green's representation formula for is 

(6.16) @\{x) G{x,Ti)C£aijp';hje'^^)dVg. 

Jm j 

Next we claim that 

\@'}+'^^M,-c'} + log{pfhi{p,))+2KY{pk,Pk)m'^\ 

(6.17) < Co(^,)4"-2-^ 

where 5 > is small positive number, 0^ is the average of 0f on M, Co(-) IS a 
positive function such that C()(5o) ^ 0° as 5o — > 0. To see this, let x = pk in (16.161 ). 



+ I y{Pk.n)Y.aijp]hje^idVg 

JB(pt.,So) j 

(6.18) +/ G{pk,^)yaijp)hJe^^dVg 



lM\B{pk.So) - ■' y ^- ^ - ^ 

Observe that 

(6.19) p^hje^'jdVgiri) = pptje^J^'^e^rfjdri = h^e^Ui]. 
Then by the definition of v*^, we obtain 

(6.20) -i- / log - T] I l^aijp'lhje'^'dVg 

JB{pt.So) j 

= -i- / (log(£,|r]i|))£a,-,.(%^.')(£.r7i)/5('^')-^.*(^*'")^r7i 

271 jB(OAejr') j ■' 
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where, by the same kind of notations used in section |3] and section 4, 

v,'''(-) = v1(m + -)-0fto,^ + £r) 

and 

(6.21) Hl\.) = ^{e,pu + ek-)e^-^''^'''''^''^-l 
By (16.7b we have 

hI\0) = pfhipu) + 0{el), VHj\Q) = £,pf V/i,-(w) + 0{el) 

and 

lll'\r]) = p'^hi{p,) + e,pfVhi{p,) ■ r] + 0{el){\ + \r]\f. 
On the other hand by Proposition 14.11 

e'-'^"') = ^>^' (1 +<I>f (t]i) + 0(q2)(i + |t7i D^-^+S). 

Hence 

(6.22) Hl'\'n)e"''^'^^ 

= {p^hM + etp^VhM ■ r] +<I>?(r])).^' (''^ + 0(£2)(i + 1171)2-+^ 
Using (16.221 ) in the evaluation of (16.201) we have 

L , log 1^^ - ^ I L ^ijPjhje^UVg 
= + ®\{pu) -Mu-c] + \og{pfhi{pt)) +Co(5o)(ef -2-^). 



Similarly 



r{Pk,Tl)Z''UPjhje'^'dV, = r{Pk,Pk)2nm^+Co{5o){e;^'-^-^) 

B(pkA)) j 



For the final term in (16.181 ) by the following crude estimate established in |[38l 

0fW = -^^^M/t + O(l), xeM\B{pu,8^), O(l)^ooif^^0 
it is easy to see 

/ G{pu,r])Y^aijp]hje^^dV, = Co(5o)0(£f-2-5)_ 

JM\B{p,,,So) j 

Thus, back to (16181 ). 

®Hpk) = + Y""' + ®iiPk) -Mu-c] + logipf hi{p,)) 

+Y{pk,Pk)27im'^ + 0{e;^'-^-^). 

Then (IgT/l) follows. 

Next for x £M\B{pk, 5o), by ( 16.161 ) and standard estimates 
(6.23) 0f(x) = 0f + 27rG(x,w)mf + Co(5o)4"-'-^ 
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Consequently by (16.231 ) and (16.171) 

(6.24) f-Q'W = e0H2^G(^,wW+£g^, xeM\B{pk,do) 



k 



where 

(6.25) \Es,\<C,{5o)ef-^+'. 

In the sequel we shall always use E§ to represent a term bounded by the right hand 
side of (EBt . 

6.1. Proof of Theorem ll.2[ In this case m < 4. Since Jj^hie®'' dVg = 1 we write 
Pi = I pfhie'^-dVg + / pfhie®-dVg = pi + pi 

Jb{pu,5q) JM\Bipk,So) 

By (I6l9l ). (lOl) 

pi = f h]e^'dr] = f h'l{e,y)e^-^y)dy. 

Let 

/i = {/ G /; lim mf = lim m'^. }. 
Using the expansion of v^^ in ( 16.151 ) we have (since m < 4) 

(6.26) pi= [ , /i?(0)^'^' W J); + o(5o)£f / e /i . 
Now for / ^ /i we have 

(6.27) pi= [ h){0)e''-^yUy + Es,, i^h 
and 

(6.28) \pi\=Es„ i^h. 

It is easy to see from (16.121 ) and (16.71 ) that the m defined by (16.121 ) is the same as 
the one in (11.61 ). 

Combining (16.261) . (16.271) . (16.281 ) and Lemma l2!2l we have 

L2n'^'" -^-'^271 271 
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Using (16.271 ). (16.281) to change from p,!^, to pf , 

-2£(mf-2)^|+<,(4)er''-'+E&. 
Using ( 16.241 ) we have 



JM\B{pt,So) 



JM\B{pk,So) hi{pk) 



Combining terms we have 

-yp^-Va-^^ 
271 V ' '2n2n 



h \ 2k Jm\b{p,,So) hi{pu) V 



(6.30) +o{5Q)e'^-' + E8,. 

We claim that for fixed k the following limit exists: 

(6.31) lim ( Sl'^- - / e^M{GM-y(P.M))_hL_dV^ 

5o^O V 271 jM\B(pt,8Q) hi{pk) V 

Indeed, write the second integral as the sum of one integral over B{pk, 5i ) \B{pk, Sq) 
and the other over M \B{pic,5i), where 5i is chosen small enough so that in 

B{pk,5i), 

(6.32) dVg = e'^dx where \ir{x) = 0{\x\^) in . 
In the local coordinate at pk, it suffices to prove 

lim ( dl'"- - I ,2nnf,(G(p,+.,p,)-y(p,.p,)) HPk + x) \ 



exists. Since 



G{pk + x,pk) = -T^-log \x\ + Y{pk+x,pk), 
271 



we use (16.321 ) and the Taylor expansions of hi and y{pk+x,pk) — Y{pk,Pk) to obtain 



e 



2nm' 



'(G(Pt+x,Pt)-7{,P,,Pt))l^.(^^^p^yi^.^p^^dy^ 



2 



x\-"''{\ + Y,c,x' + 0{\x\^))dx. 



1=1 
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where ci , C2 are some constants. Observe that the terms with ci , C2 disappear in the 
integration and < 4 — £i for some ei > for all / G /i , thus the limits in (16.311 ) 
exists. 

The following lemma states the closeness between to p if m < 4. It will be 
used to simplify the leading terms in the statements of main theorems. 

Lemma 6.2. Let ^ p eFi (p ^ Q) such that all pf — pi have the same sign, 
then 

(6.33) pf - p,- = 0(£f m'l-mi = 0{el'-^-^). 

Proof of LemmaO Let s'l = pi - pf. Then all s'l are of the same sign. By (16.301 ) 
we have 

l^{m,-2)s] = 0{er'-'). 

i 

Thus (16.331) holds and Lemma [d!2l is established. □ 

Using (16.331) in (16.301) we can rewrite the leading term as 

-yp^-Va-^^ 

271 '2n2n 
,17; V 2n Jm\b{p,,So) hi{pk) V 

Since c,, the c'^ can be replaced by c, in (16.141 ). Theorem 1 1.2l is established. □ 

Next we establish the following lemma regardless of m = 4 or not. 
Lemma 6.3. 

(6.34) V0f(O) = 2nm^Va{Pk,Pk) + 0{e',^-^-'). 
where Vi means the differentiation with respect to the first component. 

Proof of Lemma lot 

First from (16.231 ). for x € dB{pu, Sq) 

@'l{x)=@'l{x)-ff{x) 
= %\-m\ log \x-p^\^ 2%m\y(x, p^)-f\ (x) + 0(£f -^-e) . 

For the first derivative, by standard estimates 

\D\^){x)-D\-m\\og \x-p^\-r2nm\y{x,p^)-f\(x))\ < Cef-^-e^ ; = q, 1 

for X G B(/7/t,5o). On the other hand, we recall that vf(j) = 0f(£0) +21og£;(:. 
According to Proposition 14.11 

\D'{vl{y)-UKy-p,,,) - ^^(e^y) -<^^iy - p,^,)\ < Ce!{l + \y\)'-'"+'^-' 
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for Z = 0, 1 and 13^1 < IdoEj^K Using the asymptotic behavior Uj^ for \y\ 
(12.181 ) we have 

U^{y-PUk)+log{p^hi{p,)) = -dM,-m'}log\x\+c'^ + 0{ef-^) 

where \x\ = ek\y\. Thus by ( 16. 17b for / = we have 
(6.35) 

^j'{x) = 2Km^{Yix,pt)-r{Pk,Pk))-fi{^) + 0ie^-'), xeB{pk,25o)\B{pk,5o/2). 
Thus the comparison on / = 1 and (16.51) yield (16.341) . Lemma [63] is established. □ 

6.2. Proof of Theorem ll.3[ In this case m = 4, we first give a rough estimate of 
P^-Pi- 

Lemma 6.4. Let tend to Q such that all p^ — Qj have the same sign. Then 

(6.36) m'}-4 = 0{e^-'), p^ - Pi = 0{e^,-') V/ G /. 

Proof of Lemma [6l4t Since — )• 2, m = 4, then all ni^ — ^ 4. Recall that = 
27raf + 0{elr^), The rest of the proof is the same as that in the proof of Lemma 
16.21 Lemma [63] is established. □ 

By ([H 

Y,{di\ogh\{0) + MKpk))of = 0{el\oge,'), 

i 

then by ( [677] ) 

(6.37) f (^aKlog/j, + 0f)(w))af = 0(e2)ioge-i, / = 1,2. 

The distance from p* to Fi can be computed as follows: 

pf= / pfV®W„+ / p^hie^dVf,. 

■JB(pt,So) JM\B{pt,So) 

By (16.361) and Theorem 14.31 the second integral is 0{e^), this is the same as the 
computation for the single equation [,12,| . Therefore 

(6.38) pf=/ pfhe'^'dVg + Oiei) 

■JB{pt,SQ) 

= [ hU^-e^--^-d7i + 0is^) 

= I ^Hl-\r])e^-'^^^dr] + 0{el) 

= f ^phv{pt)e"-^'^^dri+ f ^ W\0M'e^-(^Uri 

JB{Q,Soet') JB{0,Soe[') 4 
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the second term on the right hand side, we first use the definition of the in 
(6.39) AHl'^O) = pfhi{p,) ( - lK{pk) + %n 



The first integral on the right hand side of the above is iKO^ + 0{s^). To evaluate 
the second term on the right hand sid( 
(leilTT) . (16361) . (lO) and (16341) to have 

hi{pk) 

+ 167r^g^-Vi7(/.,,w) + 647r2|Vi7(;.,,w)p)£f + 0(£,^-^). 

For e'^i' we use the definition of in (16.131 ) and (12.181) to have 
(6.40) pfhi{pt)e"-^'^^ =e'-\ri\-^ + 0{\ri\-^-^), |t]|>1 

for some 5o > independent of k. Using ( 16.39! ) and (16.401 ) in the computation of 
the second integral of pf^ we have 

pf = 2Ka^ + 27ie^loge,'b'le^ + 0{el). 

where 

b\ = \ (^^i^-2K{p,) + ^n + l6n^^VrY{p,,p,) + 64K'^^ 

= \[^ ^oghiipk) - 2K{pk) + 871 + \Vloghi{pk) + ^7iViY{pk,Pk)\^ 

Consequently, 

= e} log e,:' (£(4 - 2n,^)I^e^.) + 0(el). 

i 

i 

Theorem 1 1.3 1 is established. □ 

Remark 6.1. Even though there is a cut-ojf function in the definition ofx, changing 
the domain of this cut-off function does not change ^iY{Pk,Pk)- 

Proof of Theorem ll.lt 

By going through the proof of Lemma 16.31 using Theorem 11.21 for p 7^ 2 or 
Theorem 1 1.3 I f or p = Q instead of Proposition 14. 1 l one sees easily that 

0f (x) = 27im'}iYipk+x,p,) - YiPk,Pk)) -4- ff{x) + 0{ef-^) 
if p ^ Q. On the other hand 

(^>f (x) = 27im'}iYipk +x,pk) - YiPk,Pk)) - c\-ff{x) + 0{el\oge^^) 
if P = 2- Correspondingly 

(6.42) V(/)f (0) = 27rmfVi7(p,,a) + 0{ef-^) 
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if p / 2 and 

(6.43) V^^iO)=27im'^ViY{p,,p,) + Oie^loge,;') 

if P = 2- Theorem 15. 1 1 yields 

|£(5Klog/^?)(O) + 5/0f(O))af| <C4"*-2 if p / e 

i 

and 

|£(5,(log/jf)(O) + 5,0f(O))a,*^| <C£2iog£-i ifp = 2. 
(' 

Also the proofs of Theorem 1 1.21 and Theorem II .3 I give 

pf - Pi = Oief-^) = Oie'r^), m < 4 

and 

pf-p, = 0(£f*-2)logq7' =0(£flog£,-i), m = 4. 

So 

' 271 \ 0{e^loge^^) m = 4. 

By (l6!42l ) and (l643l) we obtain (fTTTT ) and (fTlIl) in Theorem O Theorem O 
is established □ 
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